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Abstract 

r~| . We present a new model of D-term dominated chaotic inflation in supergravity. The F-flat 

direction present in this model is lifted by the dominant D-term, which leads to chaotic inflation 
Q-i| and subsequent reheating. No cosmic string is formed after inflation because the U(l) gauge 

(— I ' symmetry is broken during inflation. The leptogenesis scenario via the inflaton decay in our D- 

term chaotic inflation scenario is also discussed. 

^> . PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



Among many types of inflation models proposed so far, chaotic inflation is special in 
that it can takeplace around the Planck time and make the universe large enough to avoid 
the recoUapse unless the universe is open at the beginning. On the other hand, other 
inflation models occur typically at much later times so that they suffer from the flatness 
(longevity) problem j^, namely, why the universe lives so long starting around the Planck 
scale till the low energy scale without the collapse. Furthermore, other types of inflation 
except chaotic and topological inflation also suffer from the initial condition problem [I, 0] , 
that is, why the inflaton field is homogeneous over the horizon scale and takes a value which 
leads to a successful inflation. 

For the analysis of a chaotic inflation model, one would need to consider the supergravity 
which would govern the dynamics of the early universe j^]- It, however, has been long 
considered challenging to realize chaotic inflation in supergravity simply because the F- 
term potential of a scalar field in supergravity has an exponential growth which prevents 
an inflaton from having an initial value much larger than the reduced Planck scale Mp ~ 
2.4 X 10^^ GeV. Chaotic inflation is still possible in supergravity, however, by allowing a 
non-minimal Kahler potential [l, [6|, even though it would be hard in general to justify a 
specific form of Kahler potential unless some symmetry such as the Nambu-Goldstone-like 
shift symmetry is introduced as done by Kawasaki, Yanagida and one of the present authors 
(M.Y.) for the natural chaotic inflation in supergravity [7| ^ . 

Another possibility to circumvent the above difficulty stemming from the F-term is to use 
the D-term dominated inflaton potential because a D-term does not have an aforementioned 



dangerous exponential factor. In the conventional models of D-term inflation [12|, Il3|, the 
energy density is sourced by the constant Fayet-Iliopoulos (FI) term ^ in D-term while the 
inflaton trajectory follows the F-flat direction whose slope is induced by the one-loop correc- 
tions. Since the one-loop corrections cannot exceed the tree level potential energy density, 
the whole potential energy density is at most <^ so that the inflation cannot start 
from the Planckian energy scale in such conventional D-term inflation models dominated 
by the FI term. They thus necessarily give rise to hybrid inflation. On the other hand, 
in the model of D-term chaotic inflation proposed by two of the present authors (K.K. & 
Y.M) [3], the gauge non-singlet inflaton field identified with the (almost) F-flat direction 
is automatically lifted by the dominant D-term. Such a direction can thus naturally lead 
to D-term dominated chaotic inflation. It however cannot lead to the successful reheating 
in the original model because the gauge symmetry under which the inflaton is charged is 
unbroken and the charge conservation prohibits the coupling between the inflaton and the 
standard model fields. In this paper, we present a simpler model of the D-term chaotic 
inflation in supergravity which can induce the successful reheating where the gauge symme- 
try is spontaneously broken by the non-vanishing vacuum expectation value (VEV) of the 
inflaton. As a concrete illustration of a new D-term chaotic inflation scenario, we discuss 
the leptogenesis [l^ by the inflaton decay into the right-handed Majorana neutrinos [sl, [3] 
accompanied by the sufficient reheating. 

In the next section, we present a new model of chaotic inflation in supergravity and 
discuss its dynamics and primordial fluctuations. In Sec. Illlt we discuss the leptogenesis via 



^ The application of such a shift symmetry enables us to realize not only chaotic inflation 0, S] but also its 
variants in supergravity 0, [l^ and superstring 11 1 . 
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TABLE I: The U{1) x charge assignments for the superfields. 



the inflaton decay. The final section is devoted to the summary and discussion. 



II. NEW D-TERM CHAOTIC INFLATION MODEL IN SUPERGRAVITY 

We introduce three superfields S,X,X charged under U{1) gauge symmetry and (global) 
U{1)r symmetry. The charges of the superfields are listed in Table [H which ensure our 
model is anomaly free ill . 

The general renormalizable superpotential for these fields is then given by 



W = \S{XX - fx^), 



(1) 



where we can set the constants A and fi to be real and positive for simplicity. Note that non- 
renormalizable terms S'(XX)" can appear in the superpotential. However, as shown later, 
I XX I ~ /i^ during infiation so that such higher terms are negligible as long as /U ^ Mp. 

This leads to the following scalar potential consisting of the F-term Vp and D-term Vd, 
along with the canonical Kahler potential = Xli minimal gauge 

kinetic function fabi^i) = Sab, 
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Vf 



Vf + Vd, 
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X + X*(XX-/i2 

-\xn\ 



\s\' + \s\') 

X + X*(XX 



/i 



(2) 



where g is the coupling constant of the U(l) gauge interaction and we do not introduce the 
FI term for simplicity. Here and hereafter we set the reduced Planck scale Mp to be unity 
and use the same symbols for the superfields and corresponding scalar fields unless stated 
otherwise. 

The minima of the F-term (the F-fiat condition) are given by 



XX-/i2 = o, ^ = 0, 
and the minima of the D-term (the D-fiat condition) are given by 

|X| = |X|. 

The global minima of the potential hence are given by 



(3) 



(4) 



5 = 0, X = /^e^^ X = /xe 



-ie 



(5) 



3 



where the phase 9 can be set to zero by the U(l) gauge transformation. 

One should notice that this superpotential and the corresponding scalar potential are the 
same as those of the conventional F-term hybrid inflation [l9| in which the gauge singlet 
field S plays the role of an infiaton while X and X remain zero during the inflation and then 
roll down to the global minima after the inflation. In order for the hybrid inflation to start, 
the field 5* has to be relatively large but smaller than the reduced Planck scale Mp due to 
the exponential factor in the F-term while X and X almost vanish [2ol |. 

The notable difference between the new D-term chaotic inflation and the conventional 
F-term inflation is the initial condition. In the new model, inflation occurs when |X| ^ 1 
or |X| ^ 1 with 5 ~ and XX ~ /i^ which almost satisfy the F-flat condition. When the 
universe starts around the Planck scale, the potential energy as well as the kinetic energy is 
expected to be of order the Planck energy density. This requires the almost F-flat condition 
because all the fields quickly roll down to the global minimum because of the exponential 
factor, as will be seen below. The almost F-flat direction is thus naturally realized around 
the Planck scale. The potential is consequently dominated by the D-term potential, and 
thus chaotic inflation takes place. 

We now investigate the dynamics of the new D-term chaotic inflation model in details. 
Despite the factor of F-terms, due to the presence of the relatively small but non- vanishing 
D-terms, the actual infiaton trajectory is slightly deviated from the exact F-flat direction and 
given by solving the equations (1) dV/dS* = dV/dX* = or {2)dy/dS* = dV/dX* = 0, 
depending on the initial conditions ((1) for \X\ 1 or (2) for \X\ ^ 1). Note that the 
system is invariant under the interchange of X and X so that both solutions (trajectories) 
lead to the essentially same dynamics. We therefore concentrate on the first trajectory given 
by X = X{X), S = 0, and we call this trajectory T. 

We here confirm that inflation indeed occurs along this (almost F-flat) field trajectory. 
For this purpose, we first evaluate the mass terms of the field 5* along the trajectory T, 
Vij\T(pi(pj with V^ij = d'^V/ {d(f)*d(f)j). Here represents S, X, or X. The suffix T represents 
the evaluation along the trajectory T. Then, the mass matrix of the fields S, Vijlr, is given 
by 

Vss\t ^ A^e^dXl^ + lxn, 

VsxIt = 0, ^sxIt = 0. (6) 

The effective squared mass of the field S is much larger than the Hubble parameter squared 
^ (7^|X|^/2 unless the constant A is exponentially small, which makes S quickly go to 
the zero. As a result, we can safely set S to be zero and we can discuss the dynamics of the 
infiaton based on the following potential, 

Ves{X,X) = V{X,X, 5 = 0) = A'e^lXX - /I'p + y (jXp - \X\y. (7) 

By use of the U{1) gauge symmetry, we can, for instance, make the field X real without loss 
of generality, so that the ImX rapidly goes to the zero because the effective mass squared of 
the imaginary part of X is given by 'm'^^x — A^e^X^. We therefore consider the following 
effective potential, by redefining the fields X = y/2IleX, X = a/2 Re X (we take both X 
and X to be positive for definiteness) and (x' = -\/2yU, 

VMX,X) = ^e^ {XX - + ^ - x')" (8) 
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with K = (X^ + X )/2 and the canonical kinetic terms. 

The dynamics of the inflation can be discussed based on the above potential. As men- 
tioned before, we consider the case that X ^ 1 initially, which implies X <C 1 for /i' <^ 1. 
For such a case, the trajectory (T) is characterized by the condition dV/dX = which is 
equivalent to 



.'2 „-K 



(X2 - X') X 



f X + f X (XX - 

^-i^{x'-r). (9) 

Here we have used the fact that XX — fi'"^ = 0{e~^) and X ^ X. The F-term contribution 
to the potential is then estimated as 

2X + X {XX - fi'^) V 



< / (X^ - X^) « = f {X' - X') " (10) 

for X ^ 1 and X ^ 1. The potential thus is dominated by the D-term during inflation. 

We also consider the mass matrix of the fields X and X, VijlriV^ij = d'^Vcs/ {d(f)id(f>j)), 
whose elements are given by 

y e^x' + g^ (X2 - X') X' + ^ (sX^ - x' 
^\kv^ , / f^2 —A ^ , ^2 , —2^ 



Vxx 


It ^ 


Vxx 


It ^ 


V,xx 


It ^ 



(x^-x')^(i + x2 + x')-/xx. 



2 2 



up to the order of (9((e'^)°).^ The effective squared masses, say m^, of the fields X and X 
are given as the solutions of the following equation 



- (Vxx + VxxW + VxxVj^j^ - V\-^ = 0, (12) 



where 



V,xx + I^xxIt ^ ^A^e^ [X' + X') ^ iA2e^X^ 
3 



Kxx^Xx - ^xxIt ^ 4^?'A^e^^^ (13) 

up to the order of O(e^). Here, we have used the approximation that XX — /^'^ = 0{e~^) 
and X ^ 1 ^ X. The effective squared masses are then approximately given by 

^ \x'e^X\ \g^X^ « iJ^ ^ Vz>/3, (14) 



^ More precisely, we have used the approximation that \XX — /^'^| <C -'^/-'S^, XX. 
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where if is a Hubble parameter. The inflaton field in the new D-term chaotic inflation under 
discussion corresponds to this effectively massless mode. This light squared mass vanishes 
for (yf = as expected, reflecting the exact F-flat direction. 

Since V^xx ^ ^,xx and X ^ X, the inflationary trajectory is given by the minimum of 
the fleld X represented by dV/dX = 0, and we can write the minimum of X as a function 
of X, X = X (X). The fleld trajectory governing the inflation dynamics therefore can 
be parameterized by the fleld X which we call an inflaton.^ Then, by inserting the above 
relation into the effective potential, we deflne the reduced potential Vj.{X) as 



mx) = v,s{x,x"\x)) 



9' 



X' 



(15) 



As explicitly shown in Ref. [2lJ, when there is only one massless mode and the other modes 



are massive, the generation of adiabatic density fluctuations as wells as the dynamics of the 
homogeneous mode is completely determined by the reduced potential Vt:{X). Indeed, the 
equation of motion for the homogeneous mode of the inflaton X along the rolling direction 
(T) is approximated as 



X + 3HX 



cS 



dX 



X + 3HX 



dX 



0, 



(16) 



where the dot represents time derivative. Thus, the dynamics of the inflation with the 
inflaton X can be estimated by using the reduced potential K(^) as long as the inflation 
dynamics follows the minimum of X. 

Next, we evaluate the primordial density fluctuations in the longitudinal gauge. The 
equations of motion for the perturbation SX and 5X are given by [2^ 



6X + 3H6X -6X + VxxItSX + VxxW 

(X 

5X + 3H5X -5X 



6X 



dVr 



cS 



dX 



dX 



X 



X 



4X<I>, 



(17) 



where $ is the gravitational potential. We hereafter use the same symbols X and X for both 
the homogeneous modes and the full flelds for the notational brevity unless stated otherwise. 
We are interested only in the adiabatic density fluctuations characterized by the condition 



5X _ 6X 



dX"\X) 



(18) 



where we have used X/X = dX {X)/dX. Since the relation dVes{X,X {X))/dX = 
holds for any X in the regime of our interest, we flnd 



_d_ 
dX 



^{x.JTix)) 



Vxx + y- 



dX 



XX 



dX 



0. 



(19) 



^ Note here that the effectively massless field trajectory parameterized by the inflaton field X is different 
from the X direction with the mass V xx ^ ■ 
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Then, using this relation, the equation of motion for the perturbation SX can be rewritten 



as 



V2 Y Y 

5X + 3H5X - --5X + — 
a'' 



,xx 



XX 



SX = -2 



dVr 



cS 



dX 



(20) 



Taking into account this relation and 



dX^ 



dX 



dX 



V- 



,xx 



y,xxVxx - y 



,xx 



V- 



XX 



(21) 



(22) 



the equation of motion for the perturbation 6X becomes 

d'^V dV 
6X + 3H6X - -^6X + -—^6X = -2—^$ + 4X$. 

dX ^ dX 

The perturbation bX is hence completely determined by the reduced potential V^{X^. Note 
that (P'V^I dX"^ is the effective mass squared along the rolling direction given by (9^5 / dX = 0, 
and this rolling direction actually coincides with the eigenvector of the effectively massless 
mode of A. 

On the other hand, using the adiabatic condition, the gravitational potential is described 
only by 6(f)i as 



H 



a2 



2M^ V 



( X5X - X6X 



dX'" 
dX 



(23) 



Consequently, the relation 



H 



X^ + X 
2M2 



X2 



2Ml 



djr 

dX 



(24) 



leads to the gravitational potential in the long wave limit 



dt\X J 



(25) 



This expression of the gravitational potential coincides with that of the single field inflation 
with the reduced potential Vr{X). We result calculate the density fluctuations of 

our inflation model based on the reduced potential V^(X) ~ g'^X'^/8, which implies that the 
gauge coupling g should be (7 ~ 10~^ in order to explain the primordial density fluctuations. 

After inflation, the inflaton oscillates around the global minimum X = X = fi' and S* = 0. 
Note that S remains to vanish. Though the inflaton rolls down almost along the direction 
of X during inflation, the trajectory after inflation is curved so that both X and X oscillate 
around the global minimum. The mass matrix around the minimum is given by 



Vxx\m 

y,Xx\M 

VjcxIm 



1 



2,./2 



2,./2 



-A^e^/x'^ 



2 12 
9 , 



(26) 
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where the suffix M represents the evaluation at the global minimum M. Hence, the effective 
squared masses m\ are given by = X^e^fi'"^ with the eigendirection X+ = (X + X) / \f2 
and = 1g^[i''^ with the eigendirection X_ = {X — X)j 

III. REHEATING AND LEPTOGENESIS IN NEW D-TERM CHAOTIC INFLA- 
TION 

Let us now discuss the issues of the reheating and the baryon asymmetry of the universe 
in the new D-term chaotic inflation through a concrete example, namely, the non-thermal 
leptogenesis scenario via the inflaton decay which produces the the heavy Majorana neutrinos 
TVj non-thermally [l6| . 

We consider the right handed Majorana neutrinos in addition to the Minimal Supersym- 
metric Standard Model (MSSM) fields in the superpotential 



W 



\S{XX - /i^) + XI OLiXXNiN^ + X h^^.NiLjHu + Wj 



MSSM, 



(27) 



where the subscripts i and j represent the generation indices, h^j is the Yukawa coupling, Lj 
is the lepton doublet, is the up-type Higgs and PFmssm is the superpotential of the MSSM. 
The charges for various supermultiplets are given in Table [D Note that non-renormalizable 
terms like NiNi{XX)'"- and S'(XX)" can appear in the superpotential but they are negligible 
as long as /i <C 1 because \XX\ ~ /i^ during inflation. 

Taking the canonical Kahler potential, the minimal gauge kinetic function and the van- 
ishing FI term, we can calculate the scalar potential consisting of the F-term Vp and D-term 
Vd 



V 
Vf 



Vf + Vd 

K 



X {XX - /i^) + S*W + XSX + laiXNiNi + X*W 
XSX + iaiXNiNi + XVT + ^baiXXiVi + h^^jLjHu 



N*W 



\dW 



311^1 



D 



|-(|X|2-|X|2)' + V^^SSM^ 



(28) 



where represents the MSSM field and V^ssm represents the D-term concerning the stan- 
dard model gauge group. The minimum of the F-term (the F-flat condition) is given by 

Vf = ^ XX -fi^ = & S = & Ni = k MSSM F - flat condition. (29) 

On the other hand, the minimum of the D-term (the D-flat condition) is given by 

Vd = \X\ = \X\ & MSSM D - flat condition. (30) 

The global minimum of the potential hence is given by 



5 = 0, X = /ie^^ X = ^ie 



-ie 



MSSM flat condition 



(31) 
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where the phase 6 is arbitrary. As was done in the previous section, one can show that S 
and Ni go to the zeros dynamically during the inflation so that the effective dynamics is 
described by the effective potential 

Ve^iX, X) = ^e^ {XX - /x'^)^ + , (32) 

where we have used the redefined fields X = a/2 Re X and X = ^/2IleX. Thus, the 
dynamics and primordial fluctuations of the inflation are essentially unchanged even if we 
include other fields besides X and X. 

Let us now discuss the reheating and leptogenesis in this model. Note that the spon- 
taneous breaking of the gauge symmetry due to the non-vanishing VEV of the inflaton is 
crucial because, if the inflaton VEV vanishes, the charge conservation would prohibit the in- 
flaton decay in our model. The Majorana masses of right-handed neutrinos Mj are given by 
Mi = tti {XX'^ /2 = ai/i'^/2. The inflatons X and X (or equivalently X+ and X„) can decay 
into the right handed neutrinos through the Yukawa interactions if Mj < m+/2 or m_/2. 
Let us now, for concreteness, consider the case where X ^ g and the decay of X_ to the 
right handed neutrinos is kinematically prohibited. The partial decay width of the inflaton 
to the right handed neutrinos X+ NiNi, T(^x+^n,), is then given by 



1 3 

2 / //' \ 2 



The reheating temperature Tr hence becomes 

rH^0,lvr;^~10'GeV(^)(^^j-^^5^j- (34) 

with a = af where the sum i runs for only the generations of right-handed neutrinos 

which the inflaton X+ can decay into. 

The produced decay into leptons Lj and Higgs doublets Hu through the Yukawa 
interactions via 

W = hi^N,L,H^ (35) 

where we have taken a basis where the mass matrix for Ni is diagonal. We also assume 
I(^i/)i3| > \{hi^)i2\ ^ K^i')*!! (i = 1, 2, 3). We consider only the decay of A''i assuming that 
the mass Mi is much smaller than the others, M2 and M3. The interference between the 
tree-level and the one-l oop diagrams including vertex and self-energy corrections generates 



the lepton asymmetry [15|, |23|, |2J, |25 



ei = 



T{Ni ^H^ + l) + r(iVi ^Th+~l) 



in 



hyhl 
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=2,3 
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where the effective CP violation phase is 



^cff = 7 T^^ • W'j 



m^3 I Khl 
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m,^3 here is a mass eigenvalue of the left-handed neutrino mass matrix m,^ estimated by the 
seesaw mechanism 26| as, based on our assumption that the {h^)^^ is the dominant entry 
in and M3 > Mi, 

M3 



which is consistent with the mass inferred from the Super-Kamiokande experiments [27|] for 
I (/i^)33 I ~ 10-2 and M3 ~ 10^0 GeV. 

The total decay rate of A''i, F^r-^, is given by 

r^, = r{N,^Hu + i) + r{N,^H: + l) 



Stt 

^ ^l(A.)ispM, 

— (J<M^)^(^). ,3, 

Thus, for a wide range of parameters, the decay rate Tn-^ is much larger than the decay rate 
of the infiaton T(^x+^n,) so that the produced Ni immediately decays into leptons and Higgs 
supermultiplets. 

Before estimating the lepton asymmetry produced in our model, let us evaluate the lepton 
asymmetry needed to explain the observed baryon number density. A part of the produced 
lepton asymmetry (or, exactly speaking, B — L asymmetry) is converted into the baryon 



asymmetry through the sphaleron processes, which can be estimated as [28 

ub 8 riL 



where we have assumed the standard model with two Higgs doublets and three generations. 
In order to explain the observed baryon number density, 

^ ~ (0.1- 1) X 10-1°, (41) 
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we need the lepton asymmetry, 



(0.3 - 3) X 10"'". (42) 



Now we estimate the lepton asymmetry produced through the inflaton decay. For Mi > 
10^ GeV, Ml is one hundred times larger than the reheating temperature Tr. In this case, 
the produced Ni is out of equilibrium and the ratio of the lepton number to entropy density 
can be estimated as 



— ~ - e^Br 

s 2 m+ 



10-'°5eff5, 



Ml \ / Tfi \ / m+ x-i 
109 GeVy llO^ GeVy VlOio GeV/ 



-^-~(^)-'(iL)(i^)(^)', ,43, 

where is the branching ratio of the inflaton into A^i. For M3 ~ M2 ~ m+ ~ 10'° GeV 
with a2 ~ tts = 0{1), m+ is comparable to M2 and M3 so that the decay into N2 and N3 are 
prohibited kinematically or suppressed by the phase space and hence Br ~ C>{1). Note also 
that since m_ ~ 10^ GeV for g ~ 10~^ and /i' ~ 10'^ GeV, m_ < Mi, M2, M3 so that the 
decays of X„ into all the right-handed neutrinos are prohibited kinematically, as assumed 
previously. In this case, a = ai ~ 0.1 and Tr ~ 10^ GeV (low enough reheating temperature 
to avoid the gravitino problem 2^) resulting in til/s ~ 10~'''5efr which is consistent with 



the baryon number density in the present universe. 



IV. SUMMARY AND DISCUSSION 

In this paper we have proposed a new model of D-term dominated chaotic inflation in 
supergravity. The F-fiat direction present in this model is automatically lifted by the D-term, 
which leads to chaotic inflation. The superpotential of our model was originally proposed as 
the F-term hybrid inflation where the gauge singlet field 5* plays the role of an inflaton with 
the vanishing X and X during inflation. On the other hand, we showed that another initial 
condition such that S" ~ 0, \X\ ^ 1 or \X\ 2> 1 with XX ~ /i^ can naturally occur around 
the Planck scale for a successful D-term chaotic inflation. In contrast to the previously 



proposed D-term chaotic inflation model [lj|, the inflaton can acquire the non- vanishing 



VEV after inflation which breaks the gauge symmetry spontaneously so that it can decay 
into the visible sector for the sufficient reheating'^. No cosmic string is formed after inflation 
because the U(l) gauge symmetry is broken during inflation, while such a cosmic string 
formation in the conventional D-term inflation is often problematic j3l| . 

Our model leads to the quartic potential chaotic inflation which has the tight constraints 
from the recent observations 3^. The relaxation of such constraints is possible by, for 



instance, an appropriate choice of the non-minimal gauge kinetic function such as a form 



4 



It was recently pointed out that the supergravity effects induce the inflaton decay into the visible sector 



when the inflaton acquires the non- vanishing VEV after inflation 30j , and an analogous mechanisms could 



help reheat the universe for a D-term chaotic inflation model as well. 
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/ = l + (ix|-^P + c?xl"^P {dx, d-x '■ constants) which gives a quadratic inflaton potentiaL One 
of the present authors (T.K.) also proposed the quadratic potential chaotic inflation by use of 
the FI field j33| even though a successful reheating needs more care [s^] • The consideration 
of the primordial fluctuations from a MSSM flat direction acting as a curvaton in our model 



could be another possibility for a viable quartic potential chaotic inflation model [35|, 136 
We also mention that, for a toy model using a minimal gauge kinetic function illustrated in 
our discussion, the gauge coupling g should be (? ~ 10~^ in order to explain the primordial 
density fluctuations. This value of the gauge coupling is much smaller than the standard 
gauge couplings. However, this may not be a problem because the gauge symmetry may 
be a hidden gauge symmetry, or the gauge coupling could be suppressed, for instance, by 
considering the extra dimensions. This topic would be worth further investigation. 

We have also discussed the leptogenesis scenario via the inflaton decay in this chaotic in- 
flation model, which can explain the observed baryon asymmetry for a reasonable parameter 
set consistent with the data from the neutrino experiments. 
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